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Abstract 

The presence of log-periodic structures before and after stock market crashes is considered to 
be an imprint of an intrinsic discrete scale invariance (DSI) in this complex system. The fractal 
framework of the theory leaves open the possibility of observing self-similar log-periodic structures 
at different time scales. In the present work we analyze the daily closures of four of the most 
important indices worldwide since 2000: the DAX for Germany and the NasdaqlOO, the S&P500 
and the Dow Jones for the United States. The qualitative behaviour of these different markets 
is similar during the temporal frame studied. Evidence is found for decelerating log-periodic 
oscillations of duration about two years and starting in September 2000. Moreover, a nested sub- 
structure starting in May 2002 is revealed, bringing more evidence to support the hypothesis of 
self-similar, log-periodic behavior. Ongoing log-periodic oscillations are also revealed. A Lomb 
analysis over the aforementioned periods indicates a preferential scaling factor A ~ 2. Higher order 
harmonics are also present. The spectral pattern of the data has been found to be similar to that 
of a Weierstrass-type function, used as a prototype of a log-periodic fractal function. 
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I. INTRODUCTION 



It is well known that many physical systems undergo phase transitions around specific 
critical points in the parameter space []], 0. Near these points the system is strongly 
correlated and many characteristic quantities can well be approximated by power laws, 
related to the scale-invariance of the system in that state. If we assume that <p(t) is an 
observable near a critical point, t c , as for example the susceptibility of the Ising model near 
the critical temperature, for a change of scale t — > At we have 

<j>(\t) = fjuftt), (1) 

where /i = X a since ~ t a . The power law is a solution of Eq.((TJ) for VA. 

A weaker version of the over mentioned scale invariance is the discrete- scale invariance 
(DSI) 0. In this case the system becomes self-similar only for an infinite but countable 
set of values of the parameter A. That is Eq.Q holds only for A = Ai,A2--- where in 
general A n = A n . In this case A represents a preferential scaling factor that characterizes a 
hierarchical structure in the system. The solution of Eq. (JI} can be written in a more generic 
form that accounts also for a possible discrete scale invariance: 

W)-«(gJ), (2) 

where G is an arbitrary periodic function of period 1. Using a first order Fourier expansion 
on Eq.(j2J and writing t — > \t c — t\ we obtain 

<j>(t) = A + B\t c -t\ a + C\t c -t\ a cos(u\n\t c -t\ - if), (3) 



where u = 2tt/ ln(A). The dominant power law behaviour, a hallmark of all critical phenom- 
ena, and the log-periodic corrections to the leading term are the main features of Eq. (JHJ). 

Sornette, Johansen and Bouchaud 0, ^ first pointed out how different price indices in 
the stock market show a power law increase with superimposed accelerating oscillations 
just before a crash. The remarkable fact was that the log-periodic formula ({3*j). derived for 
DSI systems, provided a very good approximation for this empirical fact. This led them 
to conjecture about the existence of a critical point in time, t c , for which the market can 
undergo a phase transition (crash). 

According to this framework the stock market is seen as a self-organized system that 
drives itself toward a critical point. Just as the Ising model has a parameter governing the 



temperature of the system and a critical temperature where the system undergoes a phase 
transition, it is postulated that the market also has such an underlying parameter which 
takes the critical value at time t c . The appearance of log-periodic oscillations has been 
related to a discrete hierarchal structure of the traders. Near the critical point, when the 
market is very compact and unstable, every perturbation can spread throughout the system: 
a common decision to sell by a certain group in the hierarchy of traders can trigger a herd 
effect, leading to a crash. This concept in a way, is similar the Sj *„ — cnticuUty 
proposed by Bak, Tang and Wiesenfeld fa], studied recently in Ref. |7|. 

Since the first paper by Sornette, Johansen and Bouchaud many physicists have been 
attracted by the idea of phase-transitions in a self-organized stock market 0, , [ll | , even 
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13, 



14J. A recent review on the subject can be found 



if criticisms have been also raised 
in Ref. Q- 

An intriguing scenario has been proposed by Drozdz and coworkers Inspired by 

theoretical consistency arguments, they found empirical evidence that short time log-periodic 
structures can be nested within log-periodic structures on a larger scale. The appearance 
of these self-similar periods, one inside the other, has been related to the underlying fractal 
nature of the DSI, giving rise to a multi-scale log-periodicity. Moreover, the existence of a 
preferential scaling ratio A ~ 2 has been pointed out for both the leading pattern and the 
related sub-structures. This last fact led them to formulate a universality hypothesis for the 
parameter A and as a consequence they fixed a priori the frequency of the oscillations to 
correspond to A = 2. In this way, the predictive power of Eq. © increases considerably [3, 



ll|]. Further evidence of embedded sub-structures has been reported recently by Sornette 
and Zhou Q,|l7|. 

Log-periodic patterns have been observed not only in bullish periods of the market but 
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3Q- 



also during the "antibubbles" , or bearish periods, that follow a market peak 
An example of these log-periodic oscillations has been documented during the long period 
of recession experienced by most of the world stock markets since the middle of 2000 [20]. 
In this period, which ends approximately (in our understanding) in the first months of 2003, 
all the most important markets world-wide are remarkably synchronized. A simple plot of 
the logarithmic indices, P{t), is provided Fig. ^ We believe this is sufficient to convince the 
reader of this behaviour. It is nothing but an expression of the growing globalization of the 
modern economy 
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FIG. 1: Time series of the logarithm of the most important indices world-wide since 2000. The 
time series have been appropriately normalized, P(t) — > ~~^jmy^~ ; where < ... > denotes the 
average over the period under consideration and a is the standard deviation. The synchronization 
of the indices during the recession period, and further on for some, is an expression of the modern 
globalized economy. The year tick in the graph marks, from now on, the position of the 1st of 
January of the year itself. 

In the present work we focus on the study of the daily closures of three of the most 
important international indices from 2000 until the end of September 2004: the DAX for 
the German market and the S&P500, the NasdaqlOO and the Dow Jones§ for the American 
market. We confirm (within our DSI model of Eq. (j2J)), the existence of a main log-periodic 
structure starting in September 2000 and ending in November 2002 for all the indices in 
Sec. |n] Moreover, for the first time, we identify a clear bearish sub-structure, starting 
around the 15th of May 2002 and ending one year later. Ongoing log-periodic oscillations 
with the same characteristic frequency, starting in January 2004, also are reported. 

We also address the question of a possible universality of the power law exponent a, 



§ The Dow Jones index has been included for historical, more that economic, reasons. In fact, it constitutes 
just a small sub-set of the S&P500 basket but, nevertheless, it is considered by many to be a good indicator 
for the health of the economy. 
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related to the trend of the time series. We then convert the S&P500, NasdaqlOO and the 
Dow Jones in Euros while the DAX is expressed in American Dollars. The conversion, while 
leaving the oscillations unaltered, seriously distorts the trends. Therefore, no universal 
characteristic can be claimed for this parameter. 

A Lomb analysis of the main structure and the relative sub-structure, presented in Sec. 11111 
reveals that the dominating frequency of both is related to a common value of A that is 
A ~ 2. A second relevant harmonic, at a frequency double that of the fundamental, is 
also present. These results confirm the fractal hypothesis of Drozdz at al. 00 

. Further 

indications pointing to self-similar log-periodicity have been found using a Lomb analysis of 



a Weierstrass-type function [22], |23(, taken as prototype of a log-periodic fractal function. 



The relevance of such a function for stock market log-periodic criticality was suggested for 
the first time in Ref. [lol ]. 



II. LOG-PERIODICITY AND SUB-STRUCTURES 

The DSI model of Eq. © is used to fit the logarithm of the DAX, NasdaqlOO S&P500 
and Dow Jones. Considering that the periodic function is arbitrary, we take the modulus 
of the cos,ne a 8 it provides a better representation of the data [11] than the cosine itself. The 
fitting procedure that we use is borrowed from Ref. According to that we express the 
parameters from A, B and C of Eq. (jBJ) as a function of a, u (or A), ip and t c by imposing 
the constraint that the cost function, % 2 , has null derivative with respect to them. Following 



the method by Drozdz at al. 
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, h|, we fix the parameter A, considering it to be a universal 
constant. In particular we choose A = 2. If this assumption turns out to be confirmed, 
then the predictive power of the model will considerably be increased. Moreover, as we are 
interested here in bearish periods only and not in predicting the most probable crash time, 
we introduce a further simplification by adjusting t c visually. In this way there are only two 
parameters left to explore: a and (p. The results of the fits are shown in Figs. El El IH and El 
for the DAX, the NasdaqlOO, the S&P500 and Dow Jones respectively. 

A log-periodic structure starting around the 1st of September 2000 and finishing in 
November 2002 is clearly evident for all the sets of data considered. Moreover a nested 
sub-structure starting around the 15th of May 2002 is also visible. Log-periodic oscillations 
also characterize the present state of the market. A possible origin of this behaviour can be 
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FIG. 2: Time series of the DAX from 1/1/2000 until 22/9/2004. A log-periodic structure of the 
approximate duration of two years and starting in September 2000 is highlighted by the solid curve 
labeled (A). In this case we fix t c =1/9/2000 in Eq. ©. A one year sub-structure is visible starting 
in May 2002, as illustrated by the dashed curve (B) (t c =16/5/2002). The dashed dotted curve 
labeled (C) (t c =26/1/2004) is related to the ongoing log-periodic oscillations. For all the fits we 
fixed A = 2. 

localized at the end of January 2004. 

At this point it is important to emphasize that the previous fits are obtained for a 
preferential scaling coefficient, A = 2. This is already a good indication of the universal 
nature of the hierarchical scaling in stock markets. 

Another interesting point regards the universality of the parameter a, related to the 
main trend of the time series. In order to have some insight into this direction we repeat 
the previous fits after the conversion of the various indices into different currencies. That is 
we transform the DAX from Euros to American dollars, and the NasdaqlOO, S&P500 and 
Dow Jones from American dollars to Euros. The results are shown in Figs. El HI and El 

It appears clear from the plots that, while the oscillatory structures are unaltered by 
the currency conversion, the trends experience a serious distortion and therefore we cannot 
extract a universal characteristic for the exponent a. On the other hand, one might wonder to 
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FIG. 3: Time series for the NasdaqlOO index from 1/1/2000 until 22/9/2004. The critical times 
and the coding are the same as in Fig. |2I The sub-structure starting in May 2002 is not as evident 
as for the S&P500 or DAX. 
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FIG. 4: Time series for the S&P500 index from 1/1/2000 until 22/9/2004. The critical times and 
coding are the same as in Fig. |2J 
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FIG. 5: Time series for the Dow Jones index from 1/1/2000 until 22/9/2004. The critical times 
and coding are the same as in Fig. [21 




FIG. 6: The DAX expressed in American dollars. Features are as described in Fig. [21 
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FIG. 7: The NasdaqlOO expressed in Euros. Features are as described in Fig. [2j 




FIG. 8: The S&P500 converted from American dollars to Euros. Features are as described in 
Fig. El 
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FIG. 9: The Dow Jones converted from American dollars to Euros. Features are as described in 
Fig. El 

what extent the market dynamics are digested in a currency other than the native currency* . 



III. A NON-PARAMETRIC APPROACH: THE LOMB ANALYSIS 

In order to justify our assumption of A = 2 we perform a non-parametric test on the 
angular frequency value of the log-periodic oscillations. Following the method proposed 
by Johansen et al. in Ref. j^jj, we analyze the time series of residuals, R(t), obtained by 
removing the leading power law trend in the logarithm of the price, P(t), according to 

= P(t)-A-BK-t\- 
y ' C\t c -t\ a v ; 

If the model of Eq. (JSJ) reproduces the behaviour of the market correctly then the residual 
dependence on in the variable ln(|t c — t\) must be a cosine function and a spectral analysis 
should reveal a high peak corresponding to the angular frequency uj. 



* The problem of a change in the currency of the S&P500 index has been addressed also in Ref. In 
this case it was argued that the main source of distortion of this index is related to the depreciation of 
the American dollar due to the feedback action of the Federal Reserve Bank. 
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FIG. 10: Lomb analysis for the DAX (solid line), NasdaqlOO (dashed line), S&P500 (dotted line) 
and Dow Jones (dashed-dotted line) during the period from 1/9/2000 to 30/12/2002. A main 
frequency around Co = 9.06 (A = 2), dashed vertical line, is clearly evident. Another important 
harmonic contribution can be seen at u ~ 2u>. In the insert the residuals, with the same coding, 
are plotted. 

Once we have obtained the residuals, shown in the inserts of Figs. [HI El and [TU1 we apply a 
spectral decomposition of these signals according to the Lomb algorithm j^. This spectral 
algorithm makes use of a series of local fits using a cosine function with a phase and provides 
some practical advantages, compared to the classical Fourier transform, when the data under 
examination are not evenly sampled, as in our case. The results of the Lomb analysis for the 
periods under consideration are presented in Figs. |HJ El and [TUl In these plots the angular 
frequency corresponding to A = 2, that is to ~ 9.06, is represented by a vertical dashed line. 

The results of the analysis show, for all the periods, a dominating peak in the vicinity of 
A = 2. This fact brings more evidence to the existence of a universal scaling factor A ~ 2. It 
is also important to underline how the Lomb analysis shows, in most of the second 
main frequency that is about two times the leading frequency. 

We can also assign a confidence level to the main peaks found in the analysis. Following 



the technique proposed in Ref. 



the ratio between the two highest peaks in the Lomb 
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FIG. 11: Lomb analysis for the DAX (solid line), NasdaqlOO (dashed line), S&P500 (dotted line) 
and Dow Jones (dashed-dotted line) during the period from 16/5/2002 to 3/5/2003. Regular high 
order harmonics are still present. As already seen from the fit, the log-periodic behaviour of the 
NasdaqlOO index is not as clear as for the other indices. The residuals are presented in the insert. 

periodogram is used to give an estimation of the significance level of the higher peak. In 
selecting the peaks for the ratio are excluded the higher-order harmonics, multiples of the 
fundamental. All the ratios found from our analysis are clustered in a range approximately 



between 3 and 6, except for the NasdaqlOO index of Fig. El Even 



or a ratio of about 3, the 



confidence level is higher than 99% assuming Gaussian noise [271 ] . If, instead, we assume 
that the noise is temporally correlated, then a fractional Brownian motion [28] with Hurst 
exponent, H, at the worst, unrealistic, case of H = 0.9, provides a confidence level which 
remains greater than 80% j27J. Hence all the peaks found by the Lomb analysis show a high 
statistical significance, with the only exception being that of the NasdaqlOO where the first 
peak has a low confidence level. 

In order to have a better understanding of the spectral patterns just found, we test 



the same method of analysis on a Weierstrass-type function 22, 23]. The Weierstrass-type 
functions are a particular solution of the discrete renormalization group equation for critical 
phenomena j^l^- Defined in the interval [0,1], these functions are characterized by a self- 
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FIG. 12: Lomb analysis for the DAX (solid line), NasdaqlOO (dashed line), S&P500 (dotted line) 
and Dow Jones (dashed-dotted line) during the period from 26/1/2004 to 22/9/2004. A main fre- 
quency along with other higher order harmonics are revealed. In this case the phase of the residuals 
for the Dow Jones differs from the other indices, reflecting a less then perfect synchronization. 

similar hierarchy of log-periodic structures accumulating at a critical point t c (t c can be 
or 1 according to particular choices of the parameters). Zhou and Sornette have shown that 
Weierstrass-type functions can provide a good approximation for the bearish period of the 
stock market starting from 2000 



1 Q 



has been used test for the Lomb 



The following Weierstrass-type function 
method: 

OO 1 

/(*) = E^exp[-A"ico S (7)]cos[A"i S in( 7 )], (5) 

where 7 G [0, 7r/2] is a parameter that fixes the oscillatory structures. If 7 = 7r/2 then the 
parameter D corresponds to the fractal dimension of the function. For 7 < 7r/2 the function 



becomes smooth and it is no longer fractal |2J, |23j but preserves the large scale log-periodic 
oscillation. Characteristic curves are illustrated in Fig. ^2 

Once we have chosen the test function, the same procedures as used for the stock market 
time series are applied to the artificial time series generated by Eq. (J5J) with A = 2. fixed. 
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FIG. 13: The Weierstr ass-type function of Eq. © for (a) 7 = tt/2, (b) 7 = 0.93tt/2, (c) 7 = 0.90vr/2 
and (d) 7 = 0. For all the plots D = 1.5 and A = 2. The sum in Eq. © has been truncated at 
N = 32 because the function does not change significantly beyond this value of N. The time axis 
has been rescaled, as the Weierstrass-type function is defined only for t € [0,1]. In the present 
plots the function is fractal only in (a). 

In this case the Lomb analysis reveals for the fractal case (7 = 7r/2), apart from the clear 
peak at the main frequency, other smaller, high frequency harmonics regularly spaced, as 
illustrated in Fig. IT21 

The higher frequency periodic peaks are a manifestation of the fractality of the function 
itself. In fact, for 7 < 1 those frequencies are absent. The high frequency harmonics, in the 
pure fractal case (7 = vr/2), are a reflection of the self-similarity of the function at different 
scales. 

The real data has similar high frequency modes that have about the same spacing as the 
ones artificially obtained with the Weierstrass-type function. We can also conjecture that 
these modes are related, as in the previous case, to self-similar structures in the time series. 
Of course, because of the lack of points and noise effects, these harmonics are not as clear 
as for the Weierstrass-type function. 

The similarity in the spectral pattern between real and artificial data is an indication 
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FIG. 14: Lomb analysis for the Weierstr ass- type function of Eq.ljKJ). The self-similarity of the 
function in the fractal case (7 = 7r/2) is reflected in the regularity of the high order harmonics. 
Once we smooth the function, the high order harmonics related to the fractality disappear while 
the dominant frequency of the log-periodic leading term is unaltered. 

of the existence of self-similar structures at different scales in stock market time series, 



ll|. The 



providing further support for the fractal framework of Drozdz and coworkers jloL 
sub-structure starting in May 2002 is a clear example of self-similarity in stock market 
dynamics. 



IV. CONCLUSIONS 



In the present work we have shown that, at least, three clear log-periodic periods have 
characterized, and still characterize, the behaviour of some of the most important indices 
worldwide since the year 2000. Moreover, one of the log-periodic structures found is embed- 
ded in a longer one, interestingly, both in the decelerating market phase. This finding sup- 



ports the hypothesis of self-similar log-periodicity proposed by Drozdz and coworkers 10, 11]. 
A non parametric analysis over these periods has also been performed. The results of the 
analysis confirm the existence of log-periodic structures. Moreover, we found further evi- 
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dence for a preferential scaling factor of A ~ 2. The presence of a higher order harmonic 
at a frequency that is double the fundamental can also be related to the fractal structure 
of the time series. A test on a Weierstrass-type function supports this hypothesis. We 
have also investigated a possible universality of the power law index a. For this purpose 
we have converted the price time series to different currencies, namely the DAX from Euros 
to American Dollars and the NasdaqlOO, the S&P500 and the Dow Jones from American 
dollars to Euros. While the log-periodic oscillations remain unaltered by this procedure, the 
trends come to be seriously distorted and no universality can be claimed. 
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